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where  cOj  is  the  driving  frequency, 

u)  is  the  natural  frequency  of  the  system, 

£  is  the  damping  ratio, 

is  the  phase  lag  between  the  driving  function  and  the 
spring-mass  system. 

To  maintain  the  overall  accuracy  of  the  mass-measuring  system  to  1 
percent,  we  will  attempt  to  keep  the  phase  lag  ^  to  a  value  less  than 
that  which  would  produce  an  error  somewhat  less  than  1  percent.  As  will 
be  shown  later,  a  value  of  phase  lag  corresponding  to  a  .6  percent  mass 
error  can  be  achieved. 

Now  consider  an  co 


and  an  u) 


and  therefore 


wi  =$ 


I  K  =  1  fit 

w  V 1 . 006M  1 . 003  \  M 


u)  j  =  1.003a;. 


Then,  substituting  in  Equation  (6), 


'J>  =  tan 


-i  2£  (1 .003)o;2  2£  (1 .003) 

^  ^  —  tan  7  — rt 


<j0 2  -  (1 ,006)0)2 


-(,006) 


For  a  damping  ratio  of  .15,  in  the  spring-mass  system,  4/  =  91. 5C 
and  therefore  the  servo  system  must  follow  to  within  1.5°, 

For  a  damping  ratio  of  .1°,  =  91.75°  and  therefore  the  servo 

system  must  follow  to  within  1.75°. 
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exactly  perpendicular  to  the  vertical,  and  that  in  fact  the  thrust 
axis  is  likely  to  wander  in  small  angles  about  the  assumed  axis. 

Also,  the  motor  with  which  we  are  concerned  in  this  paper  has  thrust 
deflectors  (jetavators)  that  produce  a  10,000-lb  side  thrust.  The 
cycling  rate  of  the  deflector  can  be  as  high  as  1  cps . 

From  the  equation  u3j  =  1.003  co ,  for  a  mass  measurement  of  0.6 
percent  accuracy,  the  frequency  must  be  determined  to  within  .3  percent. 

The  period,  T,  of  a  sine  wave  is 
(61)  T  = 


and 

(62) 

(63) 


dT  =  ~  ,.,2  1  dev  . 
cv2 


For  small  differences  the  equation  becomes 

AT  =  . 

w2 


For  Ac a  =  .003  cu,  and  co  =  102  rad/sec,  AT  becomes 


(64) 


AT 


-  2  t r  (.003) 
102 


=  190  psec 


Consider  the  system  when  the  resonant  frequency  is  102  rad/sec. 

The  force  supplied  by  the  eccentric  mass  is  then 

(65)  f(t)  =  mrco2  sin  CJt  =  (,187)(102)2  sin  102t  =  1950  sin  102t  . 

Let  the  voltage  amplitude,  E,  derived  from  the  transducer  detecting 
the  amplitude  of  the  motor’s  vertical  motion,  be  one  volt. 
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DETERMINATION  OF  ROCKET-MOTOR  MASS  BY 
MEASOREMENT  OF  THE  NATURAL  FREQUENCY 

OF  A  MASS-SPRING  SYSTEM 


By 

Benjamin  Glatt 
T est  Department 


*0* 


ABSTRACT.  A  system  is  described  whereby  the  mass  of  a 
rocket  motor  being  tested  in  a  stati:  test  stand  can 
be  determined  continuously  while  burning.  The  motor  is 
mounted  on  springs  and  driven  by  the  force  of  a  rotat¬ 
ing  eccentric  mass  so  as  to  oscillate  at  the  natural  fre¬ 
quency  of  the  mass  and  spring  system-  Given  the  spring 
constant^  the  mass  of  the  system  can  then  be  calculated 
from  the  frequency  of  oscillation.  Given  a  required 
accuracy  of  1  percent  for  the  measurement  of  the  mass 
of  the  motor;  the  specifications  for  the  control  system 
are  determined. 

A  control  system  is  designed;  dialyzed;  and  shown 
to  meet  the  requirements  of  the  system  with  the  expected 
input  signals  and  external  disturbances. 
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2.  On  page  25  of  the  subject  report,  change  Eq .  53  to  read  as  follows: 


<p  =  1_0 


-t£o>2 


and  change  the  line  of  text  that  follows  Eq.  53  to  read,  {C  where  cp  is 
the  phase  shift  and  is  the  time  constant.’’ 

3.  On  page  30,  change  paragraph  4  to  read,  ‘  ‘  To  keep  the  period  accuracy 
to  within  190  p,sec,  the  zero  crossover  points  of  the  16-cps  signal  must 
be  held  to  half  the  190  fx sec.’’  Change  Eq.  68  to  read  as  follows: 


volts 


(102) (190)10  6 
2 


volts 


.0096  volts 


M 


And  change  the  first  sentence  of  the  last  paragraph  on  page  30  to  read, 

‘ i  Therefore,  the  1-cps  signal  must  be  reduced  by  the  ratio.  .^2^  =  .0084 
or  by  -41 ,5  db.,} 

4.  On  page  53,  paragraph  4,  line  3,  delete  ‘  ‘ more  than’’;  on  line  4,  delete 
31  db’ ’  and  substitute  4 ‘41.5  db.” 
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1 


where  cl^  is  the  driving  frequency, 

Cl)  is  the  natural  frequency  of  the  system, 

£  is  the  damping  ratio, 

'F  is  the  phase  lag  between  the  driving  function  and  the 
spring-mass  system. 

To  maintain  the  overall  accuracy  of  the  mass-measuring  system  to  1 
percent,  we  will  attempt  to  keep  the  phase  lag  ^  to  a  value  less  than 
that  which  would  produce  an  error  somewhat  less  than  1  percent.  As  will 
be  shown  later,  a  value  of  phase  lag  corresponding  to  a  .6  percent  mass 
error  can  be  achieved. 

Now  consider  an  cu1 


and  an  co 


(33) 


u> 


K 


006M 


_L_  IE 

1.003  \M 


and  therefore 


(34)  cox  =  1.003CU. 

Then,  substituting  in  Equation  (6), 


(35) 


* 


tan"1  2H1 -003)0)2 

0)2  -  (1 . 006  )o)2 


-1 

tan 


24(1 .003) 
-(.006) 


For  a  damping  ratio  of  .15,  in  the  spring-mass  system,  =  91.5° 
and  therefore  the  servo  system  must  follow  to  within  1.5°. 

For  a  damping  ratio  of  .1°,  =  91.75°  and  therefore  the  servo 

system  must  follow  to  within  1.75°. 
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exactly  perpendicular  to  the  vertical,  and  that  in  fact  the  thrust 
axis  is  likely  to  wander  in  small  angles  about  the  assumed  axis. 

Also,  the  motor  with  which  we  are  concerned  in  this  paper  has  thrust 
deflectors  (jetavators)  that  produce  a  10,000-lb  side  thrust.  The 
cycling  rate  of  the  deflector  can  be  as  high  as  1  cps. 

From  the  equation  cOj  =  1.003  co,  for  a  mass  measurement  of  0.6 
percent  accuracy,  the  frequency  must  be  determined  to  within  .3  percent. 

The  period,  T,  of  a  sine  wave  is  , 

(61 )  T  = 

C 0 


and 

(62)  dT  =  dw  . 

U>2 

For  small  differences  the  equation  becomes 

(63)  AT  =  2--tAoi  . 

a.2 

For  A  Cl)  =  .003  Q)>  and  co  =  102  rad/sec,  AT  becomes 


(64) 


AT 


-  2  tt  (.003) 
■  102 


190  (j,sec 


Consider  the  system  when  the  resonant  frequency  is  102  rad/sec. 

The  force  supplied  by  the  eccentric  mass  is  then 

(65)  f(t)  =  mr  co 2  sin  cot  =  (.187)  (102)2  sin  102t  =  1950  sin  102t  . 

Let  the  voltage  amplitude,  E,  derived  from  the  transducer  detecting 
the  amplitude  of  the  motor’s'  vertical  motion,  be  one  volt. 
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(89)  Css(nT)  =  lim 

1  z  —1 


“  (.35)4  (  . 05 ) 2 


'  2 - -  (z  -  l  y 


l  + 


(20 .  2)  z  (4 .75) (z  -  .956) 


121  (z  -  l)2  (z  -  .58) 


.01  radian 


=  ,57  degrees. 

The  steady-state  error  for  the  end-of-burning  case  (system  oscillat¬ 
ing  at  16  cps)  is  then 


(90)  0  (nT)  =  lim  ~ —  i  (.031)2 

'-'SS  16  ,  z  2  (Z  -  l)3 

z  — ►  1  v  7 


1  + 


12.6  ( 3 . 1 3  )  z  (z  -  .974) 


121  (z  -  1)2  (z  -  .72) 


=  .03  radian 


=  1,7  degrees. 

It  is  seen  that  these  steady-state  errors,  when  compared  to  those 
given  on  the  bottom  of  page  15, , are  adequate  for  the  required  system 
accuracy . 


RESPONSE  TO  A  SUDDEN  MASS  CHANGE  ’ 

Given  the  open-loop  z-transfer  function  of  the  control  system 
for  the  before-burning  case 


(91) 


20(4.75  z)(z  -  .956) 

~  121  (z  -  l)2 (z  -  .58) 


the  error  function  is  given  by 

<92>  E<*>  -  r-r§b 


_ R(z) _ 

20(4.75  z)(z  -  .956) 
121  (z  -  l)2(z  -  .58) 


N0TS  CL  474  (10/63)  165 
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BACKGROUND 

The  evaluation  of  large  •Clid«4>  appellant  rocket  motors  requires 
a  knowledge  of  the  thrust  of  the  motor  throughout  its  burning  period, 
A  model  of  a  thrust  stand  if  shown  in  Tig.  1. 


S?  =.  Jfr  +  Ky  -  T  »  0 

y 

vhere  T  =  Thrust 
,  K  =  Spring  constant 

M(t)  «  Mass  as  a  function 
of  time 

y  =<  Acceleration 


Fig.  1.  Schematic  Diagram  of  Conventional 
Rocket  Test  Stand. 

Conventionally,  static  teat .stands  are  equipped  with  load  cells 
that  measure  the  displacement,  y,  and  only  the  static  information  is 
•  obtained;  i.e.,  the  equation  F  =  Ky  -  T  =  0  is  the  equation  used  to 

v 

determine  the  thrust.  However,  to  obtain  instantaneous  specific  im~ 
pulse  and  high-frequency  thrust  measurements,  the  acceleration  term 
must  be  considered. 

The  acceleration  can  be  obtained  by  placing  an  accelerometer  on 
the  rocket  motor  to  measure  the  acceleration  along  the  axis  of 
thrust.  It  then  remains  to  determine  the  mass  of  the  motor  at  any 
time  during  the  period  of  "burning.  A  method  to  measure  continuously 
the  mass  of  a  rooket  motor  vhile  burning  in  a  test  stand  is  the 
subject  of  this  paper. 


1 


MAWEPS  REPORT  jjkl 


STATEMENT  OF  PROBLEM 

Consider  a  system  whereby  a  mass  (say,  of  a  rocket  motor)  is 
mounted  on  linear  springs,  and  mounted  on  this  mass  is  a  rotating 

fl 

eccentric  mass  (Jig.  2). 


SPRING 


Fig.  2a.  Vibrating  System.  Fig.  2b.  Schematic  Diagram  of 

Vibrating  System. 

r  »  Distance  from  center  of  rotation  to  center  of  gravity 
of  the  eccentric  mass 
x  a  displacement  from  rest  position 
ax^  »  Frequency  of  rotation  of  eccentric  mass 
K  *  Spring  constant 

B  »«  Damping  factor  (inherent  in  a  practical  system) 

M  »  Mass  of  rocket  nTotor  and  rotating  motor 
m  a  Mass  of  eccentric 

The.  system  is  constrained  to  motion  in  the  vertical  direction. 
The  equation  of  motion  of  the  system  is 


(1)  M2(t)  + 


m 


X(t)  +  r  sin  to.  t 
da  x  . 


+  Bx(t)  +  Kx(t)  =»  0 


(2) 


f  2 

(M  +  m)x  -  mr  sin  co^t+Bx  +  Kx  =  0  . 
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The  Laplace  transform  of  (2)  is 


(3) 


and 


(4) 


mr 


p 

(M  +  m)  X(s)  s  +  BX(b)  s  +  KX(s)  *  — - =— 


1  + 


®L 


2 


X(s) 


1 

mr  cu^  g 


n  >  8  M  +  m  2  B  , 

1  o^2)  \  K  +  K  +  1j 


Let  oj 


K 

M  +  m  ” 


and 


2  £  B 


o>  K  ' 

then  for  the  steady  state  condition,  that  is,  as  t  ~*oo 


(5)  x(t)  a 

where 

(6) 


mr  2  2 

CD 

2  "572  772  2 '  211/2  Sin  ^®fc'  ~  ^ 


72  777772  7~27Fc 

(CD  "  Oi^  }  +  4  £  cu^  CD  |  ' 


.1  2  S  ®L  “ 

*  *  tan 


CD  -  01^ 

When  cu^  «  cd  the  system  is  being  driven  at  the  undamped  resonant  or 
natural  frequency. 

Let  o>l  a  cd  and  t  -»  00 

Then 


<7> 


x(t)  =  sin  (cut  -  |)  . 
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Equation  (j)  demonstrates  that  if  the  frequency  of  rotation  of 
the  eccentric  mass  is  the  .natural  frequency  of  the  mass -spring 
system^  the  phase  of  the  sinusoidal  motion  of  the  rocket  lags  that 
of  the  eccentric  by  90° ?  and  is  independent  of  the  damping  factor. 

Therefore t  if  this  90°  phase  shift  is  maintained  throughout  the 
test  period*  the  system  is  oscillating  at  its  natural  frequency  and 
the  mass  of  the  rocket  motor  can  be  determined  by  the  equation 

(8)  M  —  ~  m  . 

o> 

The  problem  to  be  solved  is  the  design  of  a  control  system  to 
drive  the  eccentric  mass  at  an  angular  rate  which  corresponds  to  the 
resonant  frequency  of  the  spring-mass  system. 
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SHAKER  SYSTEM 

ROCKET  MOTOR  CHARACTERISTICS 

The  method  described  in  this  paper  for  measuring  the  mass  of  a 
rocket  motor  is  applicable  to  rocket  motors  in  general $  however, 

£ 

each  system  must  be  designed  to  fit  the  characteristics  of  the  par¬ 
ticular  rocket  motor. 

The  characteristics  of  the  rocket  motor  considered  in  this 
system  ares 

1.  Weight  of  motor  before  ignition  is  25,500  lb, 

2.  Weight  of  motor  completely  burnt  is  9,650  lb. 

3.  Time  required  to  completely  expend  the  propellant  is  ap¬ 
proximately  60  seconds. 

4.  Motor  is  equipped  with  jetavators  that  deflect  the  thrust, 
thereby  producing  thrust  components  normal  to  the  main  thrust  vector 
of  approximately  10,000  lb.  The  jetavators  are  capable  of  one-cps 
operation, 

DETERMINATION  OF  SUBPENSION-SYSTEM  CONSTANTS 

The  fastening  of  the  front  of  a  rocket  motor  to  a  push  rod  in  a 
static  firing  test  stand  essentially  defines  a  mass -spring  system. 

The  designers  of  such  stands  try  to  make  the  natural  frequency  as 
high  as  possible.  Because  of  flexure  in  the  motor  cases  and  in  the 
push  rods  in  such  test  stands,  there  exists  a  limit  in  the  maximum 
natural  frequency  that  can  be  obtained.  It  has  been  determined 
empirically  that  for  a  25,000-lb  motor  the  natural  frequency  is 
between  20  and  25  cps. 
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Therefore  the  natural  frequency  of  the  system  in  the  vertical 
direction  should  be  sufficiently  below  20  cps  bo  that  the  motion 
caused  by  forces  which  can  couple  into  the  system  at  those  fre¬ 
quencies  will  be  attenuated  and  have  negligible  effect.  _lt  is  also 
desirable  to  have  a  high  natural  frequency  in  the  vertical  system, 
in  order  that  as  much  frequency  information  as  possible  can  be  ex¬ 
tracted  from  the  system.  A  natural  frequency  for  the  unburnt  motor 
of  10  cps  was  chosen  as  a  compromise  between  these  two  considera¬ 
tions  . 


SPRING  CONSTANT 

From  the  equation  to 

(9) 


if  we  Set 

2  XI 

„  CO  W 

K  =  — ~* 
6 


where 

to  »  Frequency  in  rad/sec 
K  *  Spring  constant  in  lb/ft. 


W 

M  a  —  »  liass  of  the  system  .in  slugs. 
S 


Then 


(10) 


K  .  <2  „  3.125  X  106  lb/ft  . 


NATURAL  FREQUENCY  RANGE 

The  natural  frequency  of  the  system  when  the  rocket  motor  is 
unburnt  is 


(11) 


to  a  2«  (10)  =  62.8  rad/sec  , 
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a 


and  for  the  burnt  condition 


(12) 


co 


Kg 

\IW  " 


3.125  x  10°  (32.2)  ^  / 

- - *~cj5$5' . x——~  =  102.0  rad/sec  . 


The  frequencies  in  cps  are  then  10.0  and  16.25,  respectively. 
SYSTEM  DAMPING  RATIO 


To  start  the  design  of  the  system,  the  damping  was  estimated  to 
be  15  percent  of  critical  when  the  rocket  motor  was  unburnt,  and  10 
percent  when  the  motor  was  burnt;  the  difference  in  damping  is  at¬ 
tributed  to  the  energy  absorbed  by  the  propellant  grain.  A  very  low 
damping  ratio  is  desirable  in  order  that  reasonable  vertical  dis¬ 
placement  can  be  achieved  with  little  expended  energy.  If  these  es¬ 
timates  prove  excessive  in  error,  major  redesign  is  required. 
REQUIRED  MASS  UNBALANCE 

A  limit  of  l~g  peak  acceleration  was  set  for  the  system.  Tests 
on  SNORT  (Supersonic  Naval  Ordnance  Research  Track),  had  produced  ac¬ 
celerations  of  this  magnitude  in  a  series  of  sled  tests  involving 
the  rocket  motor  being  considered  in  this  paper;  there  appeared  to 
be  no  untoward  acceleration  effects  on  the  performance  of  these 
motors  during  these  tests. 

Referring  to  Equation  (7),  and  differentiating  twice,  we  get 
the  vertical  acceleration  of  the  rocket  motor  when  excited  at  the 
natural  frequency  o>, 

k 

(13)  a(t)  a  ■-gTqgL  sin  (ojt  -  ~) 

if 

TOT  tD 

hence  the  peak  vertical  acceleration  is  ^  . 
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The  maximum  value  of  (nr)  in  order  that  the  system  not  exceed  a  peak 
acceleration  of  1  £  at  ®  =*  62.8  rad/sec  is 


2  £  Kg  2(.lg)(3.12g)(l0°)(52.2) 

10  cp.  -  '  T-  ■  (62.8)1* 

H  1.87  slug -ft. 

However,  an  existing  design  for  an  eccentric  mass  produced  a 
maximum  mr  »  .877  slug-ft,  and  was  adopted  for  use  in  the  mass  meas¬ 
uring  system. 

The  maximum  value  of  (mr)  in  order  that  the  system  not  exceed  a 
peak  acceleration  of  1  g  at  a  =*  102,0  rad/sec  is 


(15)  (nr)  ,  » 

i  5J  '16.25  cps  (io2)^ 

=»  ,186  slug -ft. 

A  device  to  reduce  the  value  of  mr  in  order  to  keep  the  maximum  ac¬ 
celeration  loading  to  below  1  g  has  been  designed  for  incorporation 
in  the  system. 

DETERMINATION  OF  DAMPING  CONSTANT  B 

The  larger  damping  factor  existing  in  the  system  before  the 
motor  propellent  is  spent  is  attributed  to  deformation  of  the  plastic 
propellent  under  vibration,  which  thus  absorbs  energy.  We  therefore 
assume  that  the  damping  changes  directly  with  the  mass  of  the  motor. 
Frequency  being  inversely  proportional  to  the  square  root  of  the 
mass,  we  can  write 

Kp 

(16)  B  »  Kj_  +  -f 
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where  and  Kg  are  to  he  determined.  We  have  the  damping  at  the 
■beginning  of  burning  and  at  the  end  of  burning.  Therefore 

(17)  14920  *  Kjl  +  {.01)1<^  , 

(18)  6120  »  !C^  +  (.00377)Kg  f 
and 


(19) 

and 

(20) 


*2 


=!  ,1 . 4l  x  10 


6  lb 


ft -sec 


Kj  *  820  — . 

The  damping  constant  at  12  end  l4  cps  is  then 


(21) 


a  12  cps 


«  820  +  — xW'~  °  10620 


(22) 


B„  ,  1  13  820  *f 

f  a  14  CpS 


1.4l  x  10 
— 195 - 


8020 


lb-sec 

ft 


DETERMINATION  OF  .MASS  UNBALANCE  FOR  VARIOUS  .VALUES 
OF  NATURAL  FREQUENCY 

The  maximum  value  of  mr  when  cu  =  14  cps  is 


(23 


<w)  .  224^  -  ,  .  S2M20  .  .382  Blug-ft, 


or  (l4  x  6.28)^  (87 J 

The  maximum  value  .of  mr  when  a)  =  12  cps  is 


(24)  <„>  .  .  SSbOoo  ,  ,797  Blug.rt. 

(12x6.28)-’  (75. **)5 

A  plot  of  the  mr  required  to  keep  the  acceleration  below  1  g,  and 
the  damping  constant  B  versus  the  frequency  of  oscillation,  is  shown 
in  Fig.  3® 
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hi » - ) - 1 - ? - - 1 - r~ 

10  li  12  !3  14  15  16 

NATURAL  FREQUENCY  (cps) 


Fig»  3.  Camping  Factor,  B,  and  Product  of  Miss  and  Center  of 
Gravity,  as  Functions  of  Natural  Frequency,  at. 
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POWER  REQUIRED  TO  DRIVE  SHAKER  SYSTEM 

The  energy  required  to  keep  a  spring -mass  system  oscillating  is 
the  same  as  the% energy  consumed  by  the  damping  In  the  system. 

The  equation  for  power  is 

(25)  p(t)  a  f(t)  v(t)  , 
hut 

(26)  f (t)  a  Bx(t)  a  Bv(t)  B  ==  damping  constant. 

Substituting,  ve  get 

(27)  p(t)  -  3v(t}2  . 

The  velocity  is  a  sine  function,  and  we  may  therefore  substitute  the 
RMS  value  of  v(t)  to  obtain  the  average  power 

(28)  P  ■  BV2  V  is  the  RMS  value  of  v(t) 

P  is  the  average  power 

The  peak  value  of  velocity  is  obtained  by  differentiating  the 
position  given  in  Equation  (7)  with  respect  to  time.  Therefore 

(29)  y  =  =  SLiiL.  •. 

2 

BIT*  03  * 

The  RMS  velocity  is  then,  — ~ — ■  ,  and  the  average  power  is 

y2  B 

(30)  P- 

|  V2  B 

The  power  required  by  the  shaker  system,  plotted  against  the 
resonant  frequencies  of  the  system,  is  shown  in  Fig.  The  values 

of  mr  and  B  corresponding  to  the  resonant  frequencies  are  taken  from 
the  curves  of  Fig.  5» 

11 
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Figure  4  shows  that  the  maximum  power  required  by  the  system 
occurs  when  mr  »  .Prrn‘  slug~ft  and.  the  peak  acceleration  is  1  g.  The 
resonant  frequency  at  maximum  power  is  11.8  cps.  The  power  is  1050 
ft~lb/sec,  The  horsepower  is  given  by; 

1050 

(51)  horsepower  <*  =  1.92 

A  5  hp  motor  was  chosen  to  drive  the  system.  The  extra  (more 
than,  double)  power  is  necessary  to  take  care  of  losses  in  the  drive 
system  and  to  bring  the  system,  up  to  speed  in  a  short  period  of  time. 
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POWER  /  FT- LB  \ 
v  SEC  / 


?ig»  4.  Power  .Required  by  Shaker  System,  Versus  Natural  Frequency 
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CONTROL  SYSTEM  REQUIREMENTS 

STEADY  STATE  ACCURACY 

A  question  that  may  be  asked  is,  Why  not  measure  the  phase  be¬ 
tween  the  driving  function  and  the  oscillation  of  the  mass -spring 
system  to  determine  the  natural  frequency  of  the  system,  and  not 
bother  driving  the  system  at  its  natural  frequency?  The  answer  to 
this  is  that  in  order  to  determine  the  natural  frequency  from  phase 
information,  the  spring  constant  and  damping  factor  must  be  known. 

The  spring  constant  can  be  accurately  determined  and  is  constant 
throughout  the  burning  of  the  motor.  However,  the  damping  of  the 
system  is  not  precisely  known  during  the  burning  period,  and  it  is 
only  at  the  natural  frequency  that  this  value  is  not  important. 

Also,  it  requires  a  great  deal  more  power  to  drive  a  system  at  a 
sufficient  amplitude  when  the  driving  frequency  is  not  near  reso¬ 
nance.  Therefore,  it  is  necessary  that  the  phase  difference  be  kept 
as  close  to  90°  as  the  specified  accuracy  of  the  system  requires. 

Tine  horizontal  accelerometer  and  its  associated  recording  equip¬ 
ment  will  be  no  more  accurate  than  1.5  percent.  In  considering  the 
evaluation  of  the  equation,  T  =  Ky  +  M(t)y,  the  value  of  M(t) 
should  be  consistent  with  this  accuracy,  and  therefore  M(t)  is  to  be 
determined  to  an  accuracy  of  approximately  1  percent.  Referring  to 
Equation  (6) 


^  a) 


ih 


if  =>  tan 
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where  m j  is  the  driving  frequency, 

t»  is  the  natural  frequency  of  the  system, 

I  is  the  damping  ratio, 

f  is  the  phase  lag  between  the  driving  function  and  the 
spring -mass  system. 

To  maintain  the  overall  accuracy  of  the  mass -measuring  system  to  1 
percent,  we  will  attempt  to  keep  the  Jihase  lag  f  to  a  value  less  than 
that  which  would  produce  an  error  of  0.5  percent. 

Row  consider  an 


(32) 

and  an  os 

(33) 

and  therefore 


ox 


/  | 

\J  1.005  M  " 


1 

1.011 


(34)  oj^  »  1.011  ox  . 

Then,  substituting  in  Equation  (6), 


A  -12{  (l.Oll)to  *„_-l  2i  (1.011) 

(J5)  +  »  to"  -  tM  — Toll 

For  a  damping  ratio  of  .15,  in  the  spring-mass  system,  f  =» 

87°55’  and  therefore  the  servo  system  must  follow  to  within  2°. 

If  the  estimation  of  damping  is  in  error  by  a  factor  of  2,  and 
|  a  .3,  then  #  «  88°57c  and  the  Bervo  system  must  follow  to  within  1°. 
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CHARACTERISTICS  OF  INPUT  SIGNAL  RESULTING  FROM  A 
SUDDEN  CHANGE  IN  MASS,  AND  REQUIRED  SYSTEM  RESPONSE 

It  is  anticipated  that  certain  failures  in  the  solid  propellant 
« 

system,  such  as  large  fragments  of  unburnt  propellent  being  blown 
through  the  nozzles,  would  cause  sudden  changes  in  the  mass  of  the 
system.  We  then  have  a  condition  of  a  step  mass  change. 

Under  normal  burning  conditions,  it  is  assumed  that  the  changes 
in  frequency  of  the  vibrating  system  are  slow,  and  that  transient 
effects  can  be  neglected.  However,  this  is  not  necessarily  true  in 
the  case  of  the  step  change.  It  Is  necessary  to  detect  these  mal¬ 
functions  when  they  occur,  and  to  know  what  effects  such  changes 
have  on  the  mass -measuring  system  as  a  whole. 

An  expression  for  the  motion  of  the  vibrating  system  when  mass 
1b  suddenly  removed  can  be  found  in  the  following  way; 

First,  assume  the  system  is  being  driven  at  resonance,  and  that 
all  transients  have  decayed  to  a  negligible  value.  At  any  time, 
t  =»  a,  evaluate  the  system  and  determine  the  position  and  velocity 
of  ■  the  mss . 

Then  consider  a  system  that  has  a  different  value  for  the  mass, 
and  is  driven  by  the  same  force  function  but  starting  at  t  a  0,  and 
whose  initial  conditions  are  the  position  and  velocity  of  the  mass 
that  was  found  for  t  =>  a- in  the  fir6t  case. 

From  Equation  (7),  the  position  .of  a  system  being  driven  at 
resonance  by  a  force,  f(t)  =•  mr  co  sin  cut,  is 

2 

.  ,,  v  mr1  o?  .  /  ,  sc\  rar  <o  .  ,  .  H\ 

x(t)  a  2~inc  Sin  “  2  **  ~~Y~  ein  "  gv  • 
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The  velocity  is  found  by  differentiating  the  above  to  give 


(36)  x(t)  =  SLiL.  cos  (cut  -  |)  . 

Now  if  the  function  is  evaluated  at  t  =>  a, 

(57)  x(a)  e  — sin  (on  -  p  - - g-  cos  cua  , 

and 

2  '  2 

;,q\  tf  \  mr  co  /  mr  co  . 

(3o)  x(a)  a  — g —  cos  (cua  -  g)  =  -g -  sin  coa  . 

Consider  now  a  sine  wave  advanced  in  time  by  "a"  second,  and 
starting  at  zero,  as  indicated  below; 


h~ *-Ja 

O 

f(t)  =  mr  co  sin  [co(t  +  a)Jp(t) 

Fig.  5*  Sine  Wave  Advanced  in  Time  by  ”aw  Second 
and  Starting  at  Zero  Time. 

The  Laplace  transform  of  f(t)  =*  mr  co^  [co(t  +  a)Jjji(t),  as  indi~ 
1 

cated  by  Aseltine  is 
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(39) 


P(s)  =  mr  o>2  C&BJ^ei n  ojt(t  ~  a) 


This  transform,  evaluated  in  Appendix  A,  is 


(40) 


mr  ay  cos  as. 


tan  as 
0) 


s  +  1 


2  2 
s  +  ay 


We  can  now  consider  the  problem  of  “the  force  f(t)  applied  to 

the  system  with  a  new  mass,  M,  and  having  the  initial  conditions 

/n\  -mr  ay 
x{0)  -  — — —  cos  as 


(1+2) 

Then 

(43) 


x  (0)  =>  + 


mr  a) 
B 


sin  as 


P(s)  -  M  X(s)s2  -  M  X(o)s  -  M  X(0)  +  B  X(s)s 
-  B  X(0)  +  K  X(s)  . 

Substituting  for  F(s), 


(44) 


[tan  coa)  . 

mr  a>  cos  os. 

-  B  +  i 

V  CO  /  J 

w  v  /  \  ^  M  BET  CO  /  \ 

=  M  X(sjs  +  — g — -  cos  (<oa)s 


X  Hh 


0) 


M  mr  o> 


sin  as  4-  B  X{s)s  +  mr  oy  cos  as  4-  K  X(s) 


and 


/tan  us  , 

mrco  cos  coa 

- - - 1  a+1 

3  0)  /  J 

sl  *f 


Murro) 

B 


cos  (cua)s 


(45)  X(s)  = 


2, 
U>  / 


Mmr  co 
B 


M  a  +  B  b  +  K 


sin  os  -  ntro>  cos  os 


s  +  B  s  4  K 
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rarco  cos  cua 
K 


tar,  aa 
oo 


M  Mb 

rr  S  •*' 

B  B 


tan  coa  +  1 


'1  8  H  2  B 

^2  j[K  s  *  k  8  +  7 


M  2  .  B  , 

7?  8  +  T?  8  +  IS 


If  a  =  6,  then 


(46)  x(s) 


aro) 

K 


M  , 

j  8  +  1 


i 


v  S  +  5  S  t  li  1  + 

a  a  i 


^  k 


0) 


2  \K 


2  A  B  , 

r?  S  +  s?  8  +  1 


K 


The  time  function  as  given  hy  Nixon  is 


(47)  (t)  « 


2  2 
*  “2 


72  2T2  7  772  2  2TI72 

(o>  -  fl>2  )  +  4  6;  ou  cu^  j  ' 


“  sin  (cot  ~  f. 


+  — sin  (oT^  1  »  t  -  to) 

*2  ^  ^ 

(Oo  V  1  -  ? 


-  OL 


2 


1  -  2  € 


M2 

B  ^72^ 

ij 


1  - 


1/2 


where 


g-Cujgt  Sin  (co  V  l"  f2  t  +  tx) 


2  4 


oo 


t,  =  tan 


co2 


2  2  * 
C0g  “  00 


-»  ’  *~r>  C0p 

*2  *  tan‘ '  s - 2 


-25  (U,2^!  -  c2 


ai  -  ajg  (l  -  2  £  ) 


2\  ' 


M  Vi  *2 

rT  OOlo  V  1  -  5 


=>  tan 


■1B^ 


,  M  » 

1  -  g-  ? 


and, 


°fe  =\/i 


If  00g  is  set  10  percent  higher  than  m  (that  is,  l.icn  =  cug)  and  4  is 
set  at  .1,  then 


0 
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(48)  x(t) 


mru>  -  . 

r=>i~7r~  4,Qco  sin 


K 


4  flcut 

49  )  +  3- 6km  0  sin  (o^t  -  232  ) 


-i  ■ 

-  5.0o»  C  sin  (uwt  +  84,3  ) 


The  exponential  terms „  when  added  vectorially  as  indicated  in  Fig. 
6*  give  a  single  exponential  term;  the  expression  for  x(t)  is  then 


(49)  x(t) 


mro) 

K 


[4.0  sin  (oat  -  49°)  +  3-5  £  ^  sin  (cc^t  +  218°) 


A  plot  of  the  two  terms  is  shown  in  Fig.  7* 


If  a  a  ~  f  then 
2o> 


(30)  x 


2  2 
cd  aVj 

2.2  ~T2  2  2]  1/2 

cu^  )  +  4  g  a)  cu^  J  ' 


—  cos  (cut 

r 


i _ e-t  v 


ojVa  -  52 


sin 


{a^VTTpt  -  t2) 


+  ~y=*x*====pj  &  sin  cu^  V  l  -  £2t 


where 


♦l  =  tan  2 - 2 


OX, 


CD 


H  -1  1  -  £2  -1  2  ^  °*2  1  “ 
=  -  tan  — qr1 - tan  -g 


^  ”~2  §  T 

ad  +  2  %  ov>  -  aip 
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SUM  OF  BOTH  CURVES 
EQUAL  TO  ZERO 
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The  time  function,  given  the  same  assumption  as  in  the  a  «  0  case, 
is 

2  f  -4  out 

(51)  x(t)  ~  p  cos  (cufc  “  49°)  +  4  2  sin  (oo^t  +  224  ) 

-tout 

+  52  sin  (o>pt) 


Adding  the  exponental  term  Vectorially  in  Fig.  8,  the  resulting  ex¬ 
pression  for  x(t)  is 


(52)  x(t) 


mro) 

K 


4  cos  (cut  -  49°)  +  3.5  2 


sin  (cu0t  -  52°) 


A  plot  of  the  two  terms  is  shown  in  Fig.  9* 

An  examination  of  Figs .  7  and  9  show  that  if  phase  measurements 
are  made  when  the  time  functions  are  zero,  the  phase  shift  ap¬ 
proaches  its  maximum  value  in  what  appears  to  be  an  exponential 
manner,  described  by  the  equation 

t 

(55)  <j>  «  1-2 

where  <|)  is  the  phase  shift  and  *s  the  time  constant. 

To  verify  the  results  of  the  analysis  of  the  system  for  a  sudden 
mass  decrease,  an  electrical  analogue  was  constructed  and  a  rapid 
change  in  inductance  was  used  to  simulate  the  mass  change. 

Given  the  following  electrical  system  (Fig.  10), 


Fig.  10.  Analogue  System. 
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the  Laplace  equation  for  the  output  voltage,,  6  (t),  is 

o 

(5*0  0n(s)  *  “o — - 5 — *  where  L,  ♦  L?  »  L  . 

°  s  LC  +  sHC  +  1 

The  Laplace  equation  for-  the  displacement  of  the  mass  in  a  mass- 
spring -damping  system  is 


Therefore  the  output  voltage  of  the  electrical  circuit  is  antilogous 
to  the  displacement  of  the  mechanical  circuit  if  we  consider  the 
inductance  analogous  to  the  mass,'  the  reciprocal  of  capacitance 
analogous  to  the  spring  constant,  and  resistance  analogous  to  the 
damping  constant.  The  values  of  inductance  and  capacitance  were 
chosen  to  produce  a  resonance  in  the  neighborhood  of  10  cps  because 
the  recording  equipment  was  suited  to  that  frequency.  Damping  of 
about  .1  to  .15  percent  critical  was  required  to  simulate  the  me¬ 
chanical  system.  With  components  that  were  available,  the  following 

system  (Fig.  11 )  was  constructed. 

*  8,5  henries 

«  120  ohms 

*  4  henries 
Rg  »  30  ohms 
0  »  32  microfarads 

Low-output 

Impedance  _ 

amplifier 


Fig.  11.  System  to  Record  Output  of  Electric 
Analogue  for  Step  Change  of  Inductance. 
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The  undamped  resonant  frequency  with  both  inductors  in  the 
circuit  is 


(56)  f  =»  - — — r  a  - - = . . -  =  :  —  a  8  CpS 

2  *  6.28  \|(12.5)02)(1Q“6) 

The  undamped  resonant  frequency  with  the  8. 5 -henry  inductor 


alone  is 


(57)  f  -  - p . - -  =  10  cps  . 

6,28  vj(8.5)(32)(lo6) 

The  damping  ratio  £  with  12. 5 -henry  Inductance  and  150  ^ 
resistance  is  given  by 


The  damping  ratio  £  with  8. 5 -henry  inductance  and  120  ^ 
resistance  is 


The  oscillator  frequency  was  adjusted  so  that  the  output  voltage 
0Q(t)  lagged  0(t)  by  90°  when  both  and  L ^  were  in  the  circuit. 

L.  was  then  removed  and  put  back  into  the  circuit  many  times  at 
random,  and  the  input  and  output  voltages  were  recorded.  Each  time 


was  added  or  subtracted,  the  phase  shifted  in  the  exponential 
manner  Indicated  by  the  analysis  and  with  a  time  constant  that 


agreed  with  the  value 


1 

f  fli. 
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REQUIRED  SYSTEM  RESPONSE 

An  overshoot  in  the  response  of  the  control  system  will  indicate 
when  a  Btep  mass  change  occurs.  But  it  is  also  desired  to  have  the 
error  settle  to  the  steady-state  value  as  quickly  as  possible.  A 
nominal  overshoot  of  10  percent  is  generally  chosen  for  a  system  in 
order  to  have  a  low  settling  time  for  a  step  input. 

It  is  therefore  desired  that  the  control  system  have  a  similar 
overshoot  when  the  exponential  transient  produced  by  a  step  mass 
change  occurs . 

UNWANTED  SIGNAL  DISCRIMINATION 

There  are  several  sources  of  signals  that  are  likely  to  be 
produced  by  the  rocket  motor  and  spring -mass  system. 

As  previously  described,  the  resonant  frequency  of  the  test 
stand  in  the  thrust  direction  is  between  20  and  25  cps .  Any  forces 
at  these  frequencies  that  couple  into  the  vertical  system  produce 
motion  that  has  an  amplitude  of  at  least  10  db  below  the  natural 
frequency  motion  for  the  same  value  of  force  amplitude.  It  is  as¬ 
sumed  that  coupling  from  the  horizontal  to  the  vertical  will  be 
small,  and  that  therefore  the  effects  of  these  high  frequencies  can 
be  neglected.  If  this  assumption  is  not  valid,  then  higher  fre¬ 
quencies  could  not  be  neglected  because  the  signals  to  the  control 
system  are  to  be  sampled,  and  aliasing  error  could  arise  even  though 
the  control  system  has  a  low,  high-frequency  cut-off. 

However,  unwanted  low -frequency  signals  must  be  attenuated. 

These  signals  arise  from  the  fact  that  the  motor  thrust  is  not 
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exactly  perpendicular  to  the  vertical,  and  that  in  fact  the  thrust 
axis  Is  likely  to  wander  in  small  angles  about  the  assumed  axis. 
Also,  the  motor  with  which  we  are  concerned  in  this  paper  has  thrust 
deflectors  (jetavators)  that  produce  a  10,000-lb  side  thrust.  The 
cycling  rate  of  the  deflector  can  be  as  high  as  1  cps . 

From  the  equation  cu^  =  1.011  ox,  for  a  mass  measurement  of  0.5 
percent  accuracy,  the  frequency  must  be  determined  to  within  1  per¬ 
cent  . 

The  period,  T,  of  a  sine  wave  is 

(61)  T  =  — 

co 

and 

(62)  dT  -  den  . 

co 

For  small  differences  the  equation  becomes 

(63)  AT  =>  Aco  . 

co 

For  Aco  =  .01  co,  and  co  =  102  rad/sec,  AT  becomes 

(64)  AT  «  =  630  psec  . 

Consider  the  system  when  the  resonant  frequency  is  102  rad/sec. 
The  force  supplied  by  the  eccentric  mass  is  then 

(65)  'f(t)  «  mrco2  sin  cot  =  (,l87)(l02)2  sin  102t  =  1950  sin  102t  . 

Let  the  voltage  amplitude,  E,  derived  from  the  transducer  de¬ 
tecting  the  amplitude  of  the  motor’s  vertical  motion,  be  one  volt. 
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Now  consider  the  force,  f^(t)  =  10,000  sin  6, 2&fc,  acting  on 
the  same  system.  From  the  frequency-response  curves  of  a  second- 
order  system  it  is  seen  that  the  response  at  1  cps  is  approximately 
13  db  below  that  at  resonance  if  the  damping  ratio  is  0.1.  There¬ 
fore  the  voltage  amplitude  derived  from  the  transducer  with  this 
force  applied  would  be 


(66) 


„  .224  x  10,000 

*T  "  1950 


1.15  volts 


where  ,224  *-13  db. 


The  slope  of  the  102  rad/sec  voltage  at  the  positive  zero 
crossover  is 


(67) 


dJE 

dt. 

t 


0 


102 


volts 

sec 


Considering  this  crossover  point,  AE  »  102  At. 

To  keep  the  period  accuracy  to  within  630  psec,  the  zero 
crossover  points  of  the  l6-cps  signal  must  be  held  to  half  the  630 
psec. 

Then 

(66)  AE  «  102  volts  *  .  •volts  *  .032  volts 

032 

Therefore,  the  1-cps  signal  must  be  reduced  by  the  ratio  * 

x  •  xp 

.028  or  by  -31  db.  The  input  to  the  control  system  will  have  to  be 
filtered  in  some  manner  to  accomplish  this  frequency  rejection. 
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CONTROL  SYSTEM 

POSITION  SERVOMECHANISM 

Consider  the  system  whereby  0 i  is  the  output  of  a  linear  posi¬ 
tion  servomechanism  driving  the  unbalanced  rotor.  Also,  consider 
the  sinusoidal  function  of  the  suspended  system  as  emanating  from  a 
rotation  that  has  a  position  0^.  Under  steady-state  conditions  the 
frequencies  of  rotation  of  the  two  systems  are  the  same.  Therefore, 
the  phase  difference  between  sin  0q  and  sin  0^  can  be  considered  as 
a  position  error,  ©Q  -  0^. 

SAMPLED  DATA  SYSTEM 

Since  the  amplitude  of  the  vibrating  system  is  not  constant, 
the  phase  of  the  two  sinusoidal  functions  can  at  best  be  compared 
every  time  the  sine  waves  cross  zero,  or  twice  every  cycle.  Com¬ 
paring  the  phase  twice  each  cycle  of  oscillation  is  equivalent  to 
sampling  the  position  error.  The  sampled-data  position  servo¬ 
mechanism  is  shown  in  Fig,  12, 


A 

G(S) 

So 

Sampler 

Fig.  12.  Sampled-Data  Position  Servomechanism. 

An  armature -controlled  dc  motor  is  chosen  to  drive  the  eccentric 

3 

mass;  the  transfer  function  of  the  motor  is  given  by 
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The  amplifier  to  drive  the  5-hp  de  motor  and  inertial  load  is  a 
60-cps  magnetic  amplifier  having  a  time  constant  comparable  to  that 
of  the  raotorj  the  amplifier  transfer  function  is  given  by 


SECOND -ORDER  SYSTEM 

A  plot  of  the  natural  frequency  of  the  system  as  a  function  of 

the  mass  is  shown  in  Fig.  13*  If  it  is  assumed  that  the  mass  change 

is  a  direct  function  of  time,  then  for  a  small  increment  of  time, 

3ct^ 

ai(t)  =  kt,  and  by  integration,  0(t)  »  -g—  ,  Therefore,  for  a  non¬ 

increasing  position  error  for  the  indicated  input,  a  double  inte¬ 
gration  iB  required  in  the  forward  loop  of  the  control  system.  The 
system  with  the  added  integrator,  magnetic  amplifier,  and  dc  motor 
is  shown  below  in  Fig.  14. 


Fig.  l4.  Uncompensated  Control  System. 

The  above  system  is  a  second-order  servomechanism,  and  without 
compensation  it  is  inherently  unstable.  Since  an  operational  ampli¬ 
fier  will  be  used  to  provide  one  of  the  integrations  in  the  system, 

4  • 

it  is  a  simple  matter  to  get  a  transfer  function  of  the  form 
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K_  (bs  +  l)(cs  +  l) 

(7!)  G(b)1  =  - - - — 

and  thereby  cancel,  one  of  the  poles  of  the  amplifier -motor  transfer 
function.  The  system  to  be  analyzed  is  shown  below  in  Fig.  15 . 


1  ' 


j  (  bs  + 1  )(cs  +  1  ) 

| 

k2  |  ; 

1 

K, 

1  ~  5  — 

i 

bs  i  1  | 

$(<35+1} 

Fig.  Ip.  Compensated  Control.  System. 

The  value  of  G ( s) ^G ( s ) a ,  as  given  by  the  manufacturer  of  the 
5 -hp  motor  and  driving  amplifier,  is 

(7o)  G(s)jnG(s)a  -  ^';o9sTT)T"i“+~iT 

Tne  mechanization  of  equation  (71)  is  shown  on  the  right  edge 
of  Fig.  17.  The  values  selected  produce  the  transfer  function 

(73)  0(3)^  .  -MlS.1  Hi*  .f.ii  . 

These  values  were  selected  so  that  the  pole  at  s  -  -10  of  G(s)a 


would  be  cancelled. 
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FILTERS  FOR  LOW-FREQUENCY  SIGNALS 

It  is  necessary  to  filter  the  input  to  the  control  system  in 

© 

order  to  provide  31  db  rejection  at  one  cps.  But  any  filter  will 
cause  phase  shift  in  the  10-16  cps  carrier t  and  the  90°  phase  re¬ 
lationship  between  the  forcing  function  of  the  eccentric  mass  and 
the  motion  of  the  rocket  motor  will  not  be  maintained.  It  is  there¬ 
fore  necessary  to  provide  the  same  phase  in  the  signal  obtained  from 
the  eccentric  mass  output.  If  identical  filters  are  used  to  filter 
the  input  signal  and  provide  phase  shift  in  the  signal  from  the  .ec¬ 
centric  output,  the  90°  phase  relationship  will  be  maintained. 

A  double  resistance -capacitance  high -pass  filter  was  selected 
for  this  purpose  because  it  is  easy  to  select  components  so  as  to 
make  both  filters  nearly  identical. 

The  cross-over  frequency  of  the  filters  was  chosen  to  be  8  cps. 
This  provides  an  attenuation  of  approximately  42  db  for  1-cps 
signals. 

SAMPLING  SYSTEM  AND  DETECTION  OF  PHASE  BETWEEN  ©in  AND  Qq 

Mounted  on  the  shaft  that  drives  the  eccentric  mass  are  two 
sinusoidal  potentiometers.  One  of  these  is  mounted  so  as  to  produce 
a  sine  wave  whose  phase  corresponds  to  the  phase  of  the  forcing  func¬ 
tion.  The  other  one  is  mounted  so  as  to  produce  a  sine  wave  that 
lags  the  forcing  fume t ion  by  90°. 

A  linear  potentiometer  is  mounted  so  as  to  detect  the  vertical 
displacement  of  the  rocket  motor. 
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The  outputs  of  these  potentiometers  are  sinusoidal  voltages. 
These  electrical  signals  are  fed  to  squaring  circuits;  that  is,  the 
sine  waves  are  amplified  and  the  top  and  bottom  of  the  wave  is 
clipped  off.  The  leading  and  trailing  edges  of  the  square  waves 
correspond  to  the  zero  voltage  value  of  the  sine  waves.  The  square 
waveforms  are  shown  in  Fig.  16. 

The  solid  square  wave  of  Fig.  l6c  has  a  phase  shift  of  less 
than  90°,  and  indicates  that  the  system  is  being  oscillated  below 
its  natural  frequency  and  that  therefore  the  driving  function  fre¬ 
quency  must  be  increased. 

The  dashed  square  wave  of  Fig.  l6c  has  a  phase  shift  greater 
than  90°,  and  indicates  that  the  system  is  oscillating  at  a  fre¬ 
quency  above  its  natural  frequency  and  that  therefore  the  driving 
frequency  must  be  reduced. 

To  obtain  the  proper  control  signal  every  time  the  sinusoidal 
motion  of  the  rocket  motor  crosses  the  zero  position,  the  following 
logic  is  employed; 

The  signal  to  increase  the  speed  of  the  control  motor  is  given 
by 

(’Jh)  =  abc+abc. 

The  signal  to  decrease  the  speed  of  the  control  motor  is  given 
by 

(75)  Sde  “  a  b  *  f  *  5  0  ' 
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Fig.  16,  Waveforms  Derived  from  Potentiometers.  Fig.  a 
in  Phaso  with  Forcing  Function;  Fig.  b  is  Lagging  Forcing 
Function  by  90  |  Fig.  c  is  Derived  from  Motion  of  Rocket 
Motor. 
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In  Pig.  16,  a,  b,  and  c  are  the  waveforms .  The  equation  above 

o 

thus  means  that  the  signal  to  decrease  the  speed  of  the  motor  occurs 
when  a  and  b  and  not  c  are  present,  or  when  not  a  and  not  b  and  c 
are  present. 

The  mechanization  of  the  logic,  and  operational  amplifier,  is 
shown  in  Fig.  17. 


SAMPLER  CONVERSION  GAIN 

1 

At  10  cps  (62.8  rad/sec)  one  radian  corresponds  to  gg-g  second. 

The  pulses  out  of  the  detection  circuit  have  amplitudes  of  25  volts. 

25 

Therefore,  a  1-radian  error  would  produce  a  vo-4k-second  pulse. 

If  a  constant  error  of  one  radian  exists,  and  if  this  error  is 

operated  on  by  one  unit  impulse  function,  a  pulse  is  produced  whose 

area  is  one  radian-second.  Since  the  sampling  circuit  produces  a 

25 

pulse  whose  area  is  jyjrg  volt -seconds  for  a  1 -radian  error,  the 
conversion  gain  Is 


(76) 


K 

(62.8) 


25 

62TF 


volt -seconds 


1  radian -6 econd 


0.4 


volts 

radian 


At  16  cps  (102  rad/sec)  a  1 -radian  error  produces  pulses  whose 
25 


areas  are  volt -seconds 


Therefore  the  conversion  gain  is 


(77) 


K  »  .245  . 

c 

(102) 
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CONTROL  SYSTEM  ANALYSIS 
ANALYSIS  AT  10  CPS  AND  AT  1 6  CPS 

The  system  described  does  not  have  parameters  that  are  constant 
in  time.  As  was  seen  In  the  previous  section,  the  gain  of  the 
system  at  10  cps  is  1.6  times  as  great  as  the  gain  at  16  cps.  Also, 
the  error  is  sampled  and  the  characteristics  of  the  sampled  system 
are  dependent  upon  the  sampling  frequency.  The  sampling  rate  of 
this  system  changes  from  20  samples  per  sec  at  the  beginning  of 
motor  burning  to  32  samples  per  sec  at  the  end  of  burning. 

Therefore  the  system  will  be  analyzed  at  the  two  extremes  to 
determine  its  performance .  It  is  to  be  observed  that  the  changes  in 
the  gain  and  the  sampling  frequency  occur  very  slowly,  compared  to 
the  time  constants  of  the  system, 
z -TRANSFORM5 

To  analyze  the  sampled-data  control  system,  the  z -transform 
will  be  used.  The  z -transform  is  defined  as 

(78)  z  =  e8T  . 

The  z -transform  describes  the  behavior  of  a  signal  at  the 
sampling  instant.  Consider  the  system  in  Fig.  18. 


/Ov  _  . 

K  (s+i) 

00 

V 

y  -  . 

s2  (s  +  (J  ) 

Fig.  18.  Sampled-Data  System  to  he  Analyzed. 
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The  z -transform  of  the  open -loop  transfer  function  can  be  fount 
by  expanding  the  function  into  a  partial  fraction  expansion,  and 
then  by  means  of  a  table  of  transform  pairs,  substituting  the  z- 
transform  for  the  s -transform  of  each  of  the  terms  of  the  expansion. 

The  expansion  of  the  open  loop  transfer  function  of  Fig.  18  is 


K  (s  +  -) 

(79)  G(s)  .  -jS.  ,  Kaa 

3  (s  *  a} 

The  z -transform  is  then^ 


”l_ 

1 

1 

1  _  1 - 

ac 

a 

c 

c  ”  a 

"2 

r  ^ 

1,1 

+  1 

_s 

8  +aJ 

(80) 


G(z)  =  Ka‘ 


f  I, 

ac 

(1  _  1' 
”  c  1 

z 

I  -  IM 

\z  -1)  £ 

E  -  eT^)J 

where  T  is  the  period  between  samples'. 

When  T  =  .05  second— that  is  when  the  sampling  rate  is  20 
samples  per  sec— and  given  ~  =  11,  G(z)  is  then  given  by 


(81) 


a(z,(20) 


Kz 


*(—  +  4.62)  +  (•—  -  4.62)] 


121  (z  -  l)d  (z  -  .58) 


,  It  can  be  seen  that,  for  —  =  1,  the  zero  comes  very  close  to 
the  pole  at  +1.  This  is  analogous  to  placing  a  zero  very  close  to 
the  origin  in  the  s -plane  when  the  effects  of  a  pole  at  the  origin 
is  to  be  offset. 

The  z -transform  when  —  =  1  is  then 

c 


(82) 


G(z) 


„  Kz  (4.75) (z  -  .95 6) 

121 7z  - 1 r  (z  -  .58) 
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When  T  =»  ,Q3  second — that  is,  when  the  sampling  rate  is  52 

1  1 

samples  per  second — and  given  —  =  11  and  —  =  1,  the  ^-transform  is 

&  c 

then 


(85)  0(z) 


(32) 


Kz  (5.13) (z 


121  (z  -  1 r  (z 

ROOT  LOCUS  ANALYSIS  IN  THE  z -PLANE 


.974) 

-  .72) 


The  imaginary  axis  of  the  s -plane  maps  into  the  unit  circle  in 
the  z -plane.  The  area  to  the  left  of  the  imaginary  axis  in  the  s- 
plane  maps  into  the  area  internal  to  the  unit  circle  in  the  z-plane. 
Therefore,  for  absolute  stability,  all  the  poles  of  the  z -transform 
must  lie  within  the  unit  circle. 

The  constant  damping  lines,  £  =  cos  0,  of  the  s -plane  map  into 
heart-shaped  curves  in  the  z-plane.  The  derivation  of  the  equation 
for  the  constant  damping  curves  in  the  z-plane  is  shown  in  Appendix 
B, 


It  is  then  possible  to  plot  the  root  locus  in  the  z-plane  and 
thereby  determine  the  gain  of  the  system  to  produce  a  characteristic 
transient  response. 

Root -locus  p3.ots  of  the  open -loop  z -transforms  given  in  Equa¬ 
tions  82  and  83  are  shown  in  Figs.  19  and  20;  the  unit  circle  and 
the  constant  damping  ratio  plots  for  £  =  0.2  and  £  =  0.4  are 
also  shown. 
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UNIT  CIRCLE 


Fig.  20.  Root  Locus  Plot  Tor  System  After  Motor  Burning 
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GAIN  DETERMINATION 

© 

To  choose  the  maximum  gain  for  the  system  to  he  consistent  with 
good  stability,  the  choice  should  be  made  for  the  case  where  the 
system  oscillates  at  20  cps.  The  gain  of  the  system  is  largest  at 
that  time,  and  because  the  system  is  being  sampled  at  the  lowest 
rate,  it  would  have  the  greatest  tendency  for  instability. 

The  open-loop  transfer  function  of  the  system  before  burning  is 
given  by  Equation  82.  The  gain  at  any  point  of  the  root  locus  can 
be  found  by^ 


(84) 


|0(*>| 


121  (z  -  l)2(z  -  .58) 

. t'."75  z(z  -  .956) 


The  values  for  the  terms  in  the  preceding  equation  are  taken 
from  the  graph  of  the  root -locus  plot  and  evaluated  for  several 
points  along  the  locus.  These  are  shown  on  the  root -locus  plot. 

Recall  that  it  was  required  to  have  an  overshoot  in  the  system 
response  to  a  step  mass  change,  and  that  the  error  should  settle  to 
the  steady -state  value  as  quickly  as  possible. 

A  damping  ratio  of  0.5  is  the  value  generally  chosen  to  meet 
the  above  requirements  when  the  input  to  a  system  is  a  step  function. 
However,  as  was  shown,  a  step  mass  change  does  not  produce  a  step 
phase  change,  but  rather  an  exponentially  rising  function.  There¬ 
fore,  a  damping  ratio  lower  than  0.5  is  chosen.  For  the  system 
before  burning,  a  value  of  |  =  .25  was  selected,  and  as  will  be 


shown  was  a  satisfactory  choice. 

The  gain  for  the  system  with  i  *  .25  is  found  by  determining 
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the  gain  at  the  point  where  the  root  locus  intersects  the  constant 
£  sa  .25  curve,  and  is 

/Qc  \  K  121  Lz^zlfiz  ~  «58)  121  (^88)  2(  .67)  pn  p  • 

(85J  \0  »  ^  •i^-(^-ri795sy-— •  =  57^  ^TTrfT  3  20,2  * 

The  gain  of  the  system  when  it  is  oscillating  at  16  cps  is 


(86) 


STEADY-STATE  ERROR 


«L6  “  Til  "  12*6 


The  steady-state  error  at  the  instants  of  sampling,  as  given  by 


7 

Jury  is 


(8?) 


^  (nT)  =  lira  -----  R(z) 
ssv  .  z  v  ' 

z-n. 


where  R(z)  is  the  z-transform  of  the  input  signal  to  the  system,  and 
G(z)  is  the  open -loop  transfer  function. 

From  Fig.  13  it  is  seen  that  the  rate  of  change  of  frequency  at 

g 

the  beginning  of  motor  burning  is  .35  rad/sec  ,  and  at  the  end  of 

2 

motor  burning  is  1.0  rad/sec  . 

The  input  to  the  system  is  then  given  by 


(88) 


e(t)  - 


kt 


which  is  9(t)  =  —  for  the  beginning -of -burning  case,  and  i8 

^2  >  v+2 

0(t)  =  g—  for  the  end -of -burning  case.  The  z-transform  for 

•  5  i  _2  zfz  +  l) 

given  by  Truxel  is  ^  kT  — - - 4-  ,  where  T  is  the  sampling  period. 

(z  -  1 V 

The  steady-state  error  for  the  beginning -of -burning  case 
(system  oscillating  at  10  cps)  is  then 
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(89)  0R  (nT)  .  lira  (.35)|  (-05)2  l 

88  10  z*l  Z  2  (a  -  l)3 


rriiogiFiEsniziiME 


=  .01  radian 


(z  -  1)  (z  -  .58)J 
»  .57  degrees. 

The  steady -state  error  for  the  end -of -burning  case  (system 
oscillating  at  16  cps)  is  then 


(90)  <2  R(nT)  6  =  liB  — ^  |  (.031)2 

88  16  z*l  Z  2  (z  -  l)5 


l^"(37f3TzT^~7yT 


1  + 


*  .026  radian 


(z  -  1)  (z  -  .72)  J 
=1.5  degrees. 

If  a  damping  ratio  of  0.2  (twice  that  estimated  for  the  end-of- 
burning  case)  is  used  in  Equation  (6),  it  is  seen  that  a  steady-state 
error  of  1.5°  is  adequate  for  the  required  system  accuracy. 

RESPONSE  TO  A  SUDDEN  MASS  CHANGE 

Given  the  open -loop  z -transfer  function  of  the  control  system 
for  the  before -burning  case 


(91) 


0(2)  .  20(^75  -?56), 


121  (z  -  l)  (z  -  .58) 
the  error  function  is  given  by 
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0 


where  R{z)  is  the  transform  of  the  input.  For  a  step  input  R(z)  = 


z  *  1 
Then 


(93) 


E(z) 


’?  -  1.58  7  +  *58  Z 


8tep  "  7  -  1.80  z2  +  1.41  z 


.58 


If  E(z)  is  expanded  in  a  power  series,  the  coefficients  of  z 
are  the  value  of  C(nT),  the  time  series.^ 

The  power  expansion  of  E(z)  for  the  step  input  is 

(94)  E(z)  =  1  + 


.22  . 

.43 

.50 

.16 

.17 

„  .22 

.07 

.08 

z 

'  2  ‘ 

3 

T~  + 

5 

►  IT  + 

7 

'  "8“ 

z 

z 

Z 

z 

z 

z 

z 

,11 

z9 

■ 

z10 

.02  , 

’  TT  + 

z 

12  + 

z 

.03 

13  ' 
z 

.01 

z 

.03 

??■ 

.02 

‘  “T5T 
z 

But  the  input  to  the  control  system  for  a  sudden  change  of  mass  is 
an  exponential  function  of  the  form 


(95)  v(t)  »  1  -  g”a4 

(96)  v(s)  a  |  -  s—T  and 


(97)  v(z)  . 

The  time  constant 

co  »  62.8  rad/sec,  £  a 
second. 

Then 


z  .  Z  /  z  \  j 

(i_~e~aT\ 

z  -  1  eaT  3  \z  -  If  1 

1  o-aTl 

>z  -  6  ' 

,  —  ,  was  shown  to  be  equal  to  .  When 
.15  and  a  ,  or  approximately  0. 


1 
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(98) 


•39  1 

1*  - 1) 

\  z  “  • 61 / 

The  error  function  with  an  exponential  input  is,  therefore,  the 

* 

•  39 


error  function  of  the  step  input  multiplied,  by 
power  expansion  is 


761 


The 


(99)  E(z) 


Exp. 


+ 

z 

2 

z 

z3 

.05 

*03 

.01 

10 

z 

11  12 
z  z 

z 


jr  *  ~ +  tt +  "9 


i?i 

c 

z ' 


The  error  functions  for  the  end -of -burning  case  are  found  in  a 
like  manner.  The  power  expansion  for  the  step  input  is  given  below. 


(100)  E(z) 


step 


z3  -  1. 

2 

,?2  z 

+  .72 

z 

z3  -  2.40 

z2  +  2.12  z 

-  .72 

,  A  .68  A 

.23 

.17 

.41 

.45 

.33 

.14 

.04 

z 

2  “ 

3  “ 

T"“ 

5  " 

7  + 

IF 

z 

z 

z 

z 

Z 

z 

Z 

.16  .20 

,  .17 

.11 

.  .04 

.ca 

.02 

0 

9  +  10 

+  11 

*  12 

73 

~  T¥ 

+  15 

z  z 

z 

z 

z 

z 

z 

z 

.03  ,  .06 

.08 

t  .08 

.06 

,  .03 

TT  *  Iff 

Z  Z 

V9 

+  20 
z 

+  21 
z 

+  22 
z 

4*  -  - 

1* 

For  the  end -of -burning  case 


is  approximately  0.1  second, 


and  time  between  samples  Is  ^  second— that  is,  T  «  .031  second. 


Therefore,  the  step  error  function  is  multiplied  by 


.27 


to 


Iz  -  .72 1 

produce  the  exponential  response.  The  exponential  power  expansion 


is 
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(101)  E(z) 


,27  ,  -38  ,35  .  .21  ,  .04 

Exp.  “  z  '  +  2  +"3  +  "TT  +  5 

*  zzzz 


.09 

z 


.16 

J 


.15 

HT 

z 


.10  .03  ,  .04  ,  .07  ^  .08  .  .07  ,  .05  ,  .03 
“9 - 10  +  “Tl+  "T2+T3  +  Tli'  +  'i5  +  'l^ 

s  z  z  z  z  z  z  z 


.  .02  .02  .03  .  .05  .  .05  .  .05  , 

■fTf,”iH+T9  +  “2o  +  "a:'>''2?+“'""  * 

z  z  z  z  z  z 

From  the  power  expansions  of  the  error  functions  the  output  can 
be  plotted.  The  output  Is  found  by  subtracting  the  error  from  the 
input.  The  transient  response  curves  for  the  four  cases  given  are 
shown  in  Fig.  21. 

From  Fig,  20  it  is  seen  that  the  overshoot  for  the  exponential 
inputs  is  about  10  percent  above  the  final  value.  The  choice  of 
damping  ratio  appears  to  be  a  good  one. 

FREQUENCY  RESPONSE 


A  knowledge  of  the  frequency  response  of  the  system  is  desir¬ 
able  in  order  that  the  effects  of  unwanted  frequency  inputs  to  the 
system  can  be  determined. 

The  value  of  G(z)  for  various  angular  frequencies  can  be  found 
from  the  z -plane  in  the  following  manner: 


(102)  •  z .  =  e8T  =  e  ^  . 

Locate  z  corresponding  to  a  frequency,  ay,  on  the  unit  circle. 

G(jdT)  can  then  be  determined  graphically.  The  factors  in  the 

numerator  of  G(jaiT)  are  the  vectors  from  the  zero  to  the  value 
1  (nT 

z  =1  2  k;  the  factors  in  the  denominator  are  the  vectors  from  the 

IcuT 

poles  to  the  value  z  =  0 u 
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O.l  0.2  0.3  0.4 


RESPONSE  TO  STEP 
CHANGE- BEFORE 
BURNING 


SEC 


RESPONSE.  TO 
EXPONENTIAL  FUNCTION- 
BEFORE  BURNING 


0.4  SEC 


RESPONSE  TO 

exponential  function- 

end  OF  BURNING 


0.4  SEC 


Fig.  21.  Transient  Response  of  Control  System 


FREQUENCY  RESPONSE 
8EF0RE  BURNING 


Fig.  22.  Frequency  Response  of  Control  System. 
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The  ratio  of  the  output  to  the  input  is  given  by 


The  expression  1  +  0(z)  can  be  determined  by  adding  one  to  the 
9o 

vector  G(z).  3—  is  then  found  by  the  division  of  vector  quantities* 

win 

The  frequency  response  found  in  this  manner  for  the  beginning- 
of -burning  case  and  for  the  end -of -burning  case  is  shown  in  Fig.  22. 

An  examination  of  the  frequency  response  shows  that  the  gain  in 
both  cases  at  1  cps  is  very  slightly  greater  than  unity,  and  there¬ 
fore  the  double  rc  filter  with  an  attenuation  of  42  db  is  more  than 
adequate  in  providing  the  31  db  attenuation  required  for  the  speci¬ 
fied  system  accuracy. 

It  should  also  be  pointed  out  that  if  an  unwanted  signal  of 
3  cps  is  present,  the  rc  filter  provides  only  23  db  attenuation, 
while  the  system  has  a  gain  of  7  db  (2.25  voltage  gain)  to  give  a 
net  attenuation  of  only  16  db. 

REACTION  TORQUE 

The  shaker  system  vibrating  in  a  vertical  direction  exerts,  a 
reaction  torque  on  the  eccentric  mass  of  the  shaker  motor. 


Fig.  23.  Rotating  System  with  the  Center  of  Rotation 
Moving  Sinusoidally  in  the  Vertical  Direction. 
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the  vertical  osoillation  when  the  system  is  vibrating  at  the  natural 
frequency  is  given  by  Equation  (7)# 

x(t)  =  sin  (cut  -  |)  . 

The  force  on  the  mass  is 

(104)  f(t)  =■  mx  . 

The  reaction  torque  is  then 

/  v2  4 

(105)  Tp  =  f(t)  r  cos  cut  =  sin  (cut  -  g-)  cos  cut 

,2  4 


cu 


1  +  cos  2  cut 


2  %  K  1  2 

In  a  manner  similar  to  that  shown  in  Brown  and  Campbell,"'  the 
transfer  function  of  the  dc  motor  will  be  derived  with  the  reaction 
torque  T^  as  part  of  the  load.  The  schematic  of  the  motor  and  load 
is  shown  in  Fig.  24. 

0(t)  =  applied-  voltage 


Re 

T, 


R 


=  armature  resistance 

=  reaction  torque  load 

J  «  armature  and  -load 
moment  of  inertia 

0(t)  =  output  shaft  position 

0,  a  counter  emf 


Fig.  24.  Schematic  Diagram  of  dc  Motor  and  Load. 

From  the  equivalent  circuit  and  the "equalities  given,  we  may 
write  the  following  equations: 

(106)  ,E(s)  a  Ra  Ia(s)  +  E^s) 

(107)  \(b)  «  K,,  s  0(b) 
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(108)  T(b)  -  Kg  Ia(a)  s 

(109)  T(s)  *  Jb2  0(b)  +  Tr(b)  . 


By  substitution,  the  following  equation  is  found; 


(no) 


Q  6 

Tr(s)  »  0(s)  s  s  + 


The  above  equation  shows  that  the  reaction  torque  can  .be  con¬ 


sidered  as  an  input  voltage  subtracted  from  the  normal  input  to  the 
motor  as  indicated  below  (Fig.  25)? 


Fig.  25.  Equivalent  System  of  dc  Motor  with 
Reaction  Torque . 

The  complete  control  system  is  shown-  in  Fig.  26. 


Fig.  26.  Control  System  with  Reactive  Torque  Input'. 

For  ©in(3)  =  0*  the  z -transform  of  the  output  for  the  system  in 
Fig.  26  is 
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RG2(z) 

Vz)  3  ”  r T g3o2<zJ 


where 


GjiCs) 


KgK5(c  s  +  l) 


(115) 


K. 

0o(b)  =»  ■  — i — . — ■ 

2W  s  (a  a  h 


(114) 


A  1 

R(b)  =  —•  TR(s) 


(115) 


R(t)  -  ^  TR(t)  - 


R  (rar)^  o)^  „  0  , 

a  '  /I  +  cos  2  cot 


{•116 ) 


R&  (mr)^  co* 


Then 


(117) 


RG2(s)  + 


s  (a  8 


and  the  time  function  as. given  by  Nixon  is 

~IJL  \  I  2  -  1 

(118)  i»p(t)  «  Rj^  Jea  +  |  - 1 +  [  ■--*  e  a 

\  /  4  to  \1  +  a  4to 


2o>  sin  ‘(2cot  -  y) 

/,  2  .  27172 

t(l  +  a  4u>  )  '  j 


The  steady  state  value  is 
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'  (H9) 


j  2  to  sin  (2tot  -  f)1 
'■to2  (l  +  4a2 n2)1^2] 


=  BK^ 


sin  (2yjt  -  y) 

2m  (i  +  4a2o)2 ) 1  /2 . 


where 


ijr  =  tan  ^  2  a  o>  . 
Sin  (2o)t  ~  f)  can  be  written  as 


( 120)  sin  (2ojt  -  t)  = 


sin  2ast 


2axo 


vV5  * 


cos  2o)t 


a  co  +  1 


A?  2 


a>~  +  1 


Then 


(121)  P0g(t).S8  -  BK^ 


and  the  z -transform  from  Jury  is 


,  sin  2mt  «  a  cos  2cot  ] 

u'r^i,'ir 77377 

7 


(122)  RG2(z)ss  -  BK^ 


T  z 


z  sin  2u>t 


(z  -  l)2  (z2  -  2z  cos  2o/T+l)(l  +  4a2o>2} 


2m 


But 

(123) 

Therefore, 


az  (z  -  cos  2mt) 


p  P 

(z^  -  2z  cos  2oyr+l)(l  +  4a  m  ) J 


2f 


sc 

«>£*> 

CD 


(124)  RGgU)^  =  BK^ 


az  (z  -  l) 


(z  ~  l)1-  (z2  -  2z  +  1)(1  +  4a2m2) 


biclI 


T  z 


(z  ir 


_ az  _ ] 

1)(1  +  4a2  to2)} 
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It  is  to  tie  noticed  that  in  sampling  a  sine  function  at  the 
frequency  of  the  function,  the  result  appears  as  a  dc  value.  The 
function  — is  the  transform  of  a  step  function. 

Z  *  -l 

For  the  before -burning  case 


(125) 


G^z)  -  G(z) 


K  z  (4.75) (z  -  .956) 
121  (z  -  l)2(z  -  .58) 


(126) 


RG2(z) 

9o(z)  =  1  +  "G^OgCz)' 


The  “steady-state  error  of  ©  (z)  from  the  final  value  theorem  is 

o 


(127)  0  (z)  =  lira  £-I_i  rg2(z) 


o'  ss  .  z 
z-»l 


1  + 


BKl 


T  z 


lim 

Z-*l 


az 


z  -  1  ~LL(z  -  l)2  (z  -  1)(1  +  4a2co2)j 


x  ~ . Kz . (4.75)(z  -  .9581 

121  (z  -  l)2(z  -  .58) 


=  0  , 

The  after-burning  case,  having  a  transfer  function  of  the  same 


form  as  the  before -burning  case,  produces  the  same  results. 


* 
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CONCLUSION 

The  system  analysis  indicates  that  the  objectives  of  the  system 
can  be  met  with  the  design  described.  The  control  system’s  steady- 
state  error  was  determined  to  be  *57°  for  the  beginning -of -burning 
case  and  1.5°  fWP  th@  end -of -burning  case }  it  was  shown  that  if  the 
damping  of  the  system  were  twice  that  which  was  estimated,  these 
steady -state  error  values  would  produo®  a  mass  error  of  less  than 
,5  percent. 

For  a  sudden  change  in  mass  the  system  responds  with  an  over¬ 
shoot  of  approximately  10  percent,  providing  an  indication  of  the 
change,  and  accomplishing  the  change  so  that  the  new  steady-state 
value  is  approached  in  a  short  period  of .time. 

The  frequency  response  of  the  system,  in  conjunction  with 
filters  in  the  feedback  loop  and  in  the  input  to  the  control  system, 
provides  adequate  discrimination  against  unwanted  low  frequencies, 
which  amplitude -modulate  the  10-16  cps  carrier.  Assumptions  are 
made  that  other  frequencies  not  generated  by  the  shaker  will  not 
have  amplitudes  significant  enough  to  affect  the  accuracy  and. 
operation  of  the  system.  These  assumptions  should  be  verified  ex¬ 
perimentally;  if  they  should  prove  to  be  wrong,  steps  would  be 
required  to  eliminate  or  discriminate  against  these  frequencies , 

The  usefulness  of  the  mass  determining  system  is  contingent  upon  the 
amplitudes  of  these  frequencies  being  insignificant. 

The  effect  of  the  vertical  oscillation  feeding  back  into  the 
control  system  via  the  eccentric  load  on  the  dc  motor  was  shown  to 
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be  zero,  and  therefore  it  would  offer  no  problem  to  the  system. 

The  analytic  study  demonstrates  the  feasibility  of  measuring 
the  mass  of  a  burning  rocket  motor  by  determining  the  natural  fre¬ 
quency  of  a  spring -mass  system. 

An  experimental  system  is  being  constructed  at  the  U.  S.  Naval 
Ordnance  Test  Station.  As  soon  as  this  laboratory  system  Is 
completed,  experiments  will  be  conducted  to  verify  the  results 
^obtained  in  this  paper.  The  proof  of  the  scheme,  however,  can 

o 

come  only  when  it  is  applied  to  a  test  stand  containing  a  burning 
rocket  motor. 
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APPENDIX  A 


Derivation  of  the  Laplace  transform  of  a  function  that  is  zero 
before  time  zero,  and  is  a  sine  function  delayed  by  a  time  **8^  for 
all  time  greater  than  zero. 

The  function  'described  in  the  title  is  written  as 
(128)  f(t)  =  sin  o)(t  +  a)  p(t)  , 

and  is  pictured  as 


f(t)  =  nr  «)2  sin  j<u(t  +  a)|(i{t) 

Fig.  27.  '  Sine  Wave  Advanced  in  Time  by  ”an  Second 
and  Starting  at  Zero  Time, 

From  Aseltine^  then  we  have  the  identity 

f  00 

(129)  <)£^jsin  o>(t-ha)  ,u(t)]  =  gas j  sin  cot  u(t-a)  <2~Sj  dt 
Therefore, 

(150}  F(s)  =  Cas  I  sin  cut  u(t~a)  <£"st  dt 
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~rCX> 

3  6aS  I  sin  cot  p(t-a)  0-S  dt  . 


Writing  sin  cut  in  the  exponential  form,  we  have 


-00 


(131) 

Integrating  gives 


x  0aS  f  Lieut  --lcot]  -st  ,, 
f(s)  a  21"]  e  -  0  e  At  . 


(132) 


F(s) 


>as 


2i 


gicut  -  st  -iait  -  st 


,oo 


10)  -  s 


-ico  -s 


Ja 


e 


as 


21 


f  ion  -  sa  --icoa  -  sa 
K  K 


ico  -  s 


-ioo  -s 


Multiply  the  numerator  and  denominator  of-  both  terms  by  the  con¬ 
jugates  of  the  denominators  -to  give 


(133) 


F(b) 


~sa  ^sa 
0  0 

2i(cu2  +  s2)  1 


f.  JLasa  •  „~icua  JLtua  r,~ 
i 0)0  +  i o)0  +  s6  -  s0 


Putting  th'e  exponential  terms  into  the  trigonometric  form  gives 

1 


(13^) 


F(s) 


2  2 
0)  +  s 


o)  cos  cua  -f  s  sin  cm 


factoring  cu  cos  tua  gives 
(135)  F(s ) 


co  cos  cm  ft 

3  tan  ob  ,  i \ 

2  2  » 

u  . ..  -t*  1 1 

.  03  1 

* 

0)  +  s 

Iob 


62 


NAVWEPS  REPORT  7?4l 


APPENDIX  B 


Mapping  the  Damping  Ratio,  ^ ,  into  the  z-Plane, 
The  z -transform  is  defined  as 


(156)  z  =  esT  . 

Let 

(137)  z  = 

and 

(138)  s  =  re1®  . 

Putting  Equation  (136)  into  the  logarithmic  form  gives 


(139)  log  R  +  itjr  a  Tre19  . 

Putting  the  term  on  the  left  side  of  Equation  (4)  into  the  polar 
form  gives 

(140)  e1  tan"1  E**  .  . 

Therefore 

(141)  ■  • 

Now 

(142)  .  4  =  -  cos'  9  , 

Therefore 

-  £ 

(143)  tan  9  =  - . .  . 

VT772 ' 

_  i  ... 

Substituting  - - -  for  tan  0  in  Equation  (141)  gives 

1-4 


(144) 


log  R  =  - l_i_ 

O 
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Tailing  the  anti-log  gives 

(145)  R  =  e 


Equation  (l45)  defines  a  curve  in  the  z-plane  for  fixed  values  of  £. 
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